Ferromagnetic polarons in Lao.sCao.sMnOs and La .33Ca .67MnO3 



G. Zheng and C.H. Patterson 
Department of Physics and Centre for Scientific Computation, 
University of Dublin, Trinity College, Dublin 2, Ireland. 
(February 2, 2008) 



Unrestricted Hartree-Fock calculations on Lai-zCa^MnOs 
(x = 0.5 and x = 0.67) in the full magnetic unit cell show 
that the magnetic ground states of these compounds consist 
of 'ferromagnetic molecules' or polarons ordered in herring- 
bone patterns. Each polaron consists of either two or three 
Mn ions separated by O - ions with a magnetic moment op- 
posed to those of the Mn ions. Ferromagnetic coupling within 
the polarons is strong while coupling between them is rela- 
tively weak. Magnetic moments on the Mn ions range be- 
tween 3.8 and 3.9 /j,b in the x = 0.5 compound and moments 
on the O - ions are -0.7 \xb- Each polaron has a net magnetic 
moment of 7.0 hb, in good agreement with recently reported 
magnetisation measurements from electron microscopy. The 
polaronic nature of the electronic structure reported here is 
obviously related to the Zener polaron model recently pro- 
posed for Pro.6oCao.4oMn03 on the basis of neutron scattering 
data. 
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I. INTRODUCTION 

The current paradigm for the electronic structure of 
manganites (Ai-^B^MnOs) with x > 1/2 is a lattice of 
Mn sites with d 3 or d 4 orbital occupancy, with the pro- 
portion of each decided by the value of x and orbital or- 
dering (OO) of the occupied e g orbital on d 4 sites. The 
corresponding double-exchange model Hamiltonian has 
been studied extensively 1-3 . Manganites with x > 1/2 
usually exhibit a phase transition 4-9 which has been as- 
sumed to be charge ordering (CO) of the d 3 and d 4 sites 
at a temperature well below the paramagnetic transition 
temperature. The conventional picture of the orbital and 
charge ordered phase with the CE magnetic structure 10 
is shown schematically in Fig. la. 

The validity of this picture for the manganites has been 
questioned recently 11 ' 12 ; its contradictions were pointed 
out by Goodenough in 1955 13 . An alternative picture, 
the Zener polaron 14 , which challenges the ideas of double- 
exchange and CO in the manganites, was recently pro- 
posed to account for the structure of Pro.6oCao.4oMn03 
determined by neutron scattering 15 . In the conventional 
CO and OO picture, ordering of the e g electron on d 4 
Mn sites is expected to induce a Jahn- Teller (JT) dis- 
tortion; the long Mn-0 bond distance in LaMnOa ex- 
ceeds 2.15 A 10,16 while the Mn-0 distance in cubic per- 
ovskite CaMnOs is 1.88 A 10 . The lesser modulations of 
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FIG. 1. Schematic illustrations of the magnetic unit cell for 
Ao.sBo.sMnOa compounds with a CE magnetic structure: (a) 
the conventional double exchange picture; (b) the polaronic 
picture; (c) a possible change in orbital order in a small fer- 
romagnetic domain when an electron is added at the oxygen 
ion site indicated by an unshaded circle. Spins are indicated 
by arrows, Mn d 4 sites by double-lobed orbitals, Mn d 3 sites 
by small shaded circles, O 2- ions by larger shaded circles and 
0~ ions by black circles. 



bond length in doped manganites with x ~ 1/2 (1.92 
and 2.06 A in Lao. 5 Cao.5Mn0 3 5 and 1.98 and 2.01 A in 
P r o.6oCao.4oMn03 15 ) suggest an intermediate valence. 
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In this Communication we 

report results of Unrestricted Hartree-Fock (UHF) cal- 
culations on Lai-^Ca^MnOs for x = 1/2 and x = 2/3. 
The UHF electronic structure for these compounds is in- 
terpreted in terms of ferromagnetic polarons containing 
two (x = 1/2) or three (x = 2/3) Mn ions. The main dif- 
ference between results from UHF calculations and the 
conventional double exchange picture is that all Mn ions 
are essentially d 4 . Consequently electrons must transfer 
from oxygen ions to every other Mn ion in the x = 1/2 
compound, resulting in an ordered array of 0~ ions be- 
tween pairs of Mn ions. The structures of the polaron 
phases for x = 1/2 and x = 2/3 are shown schematically 
in Figs, lb and 2. 



FIG. 2. Schematic illustration of the magnetic unit cell 
for Lao.33Cao.67Mn03 with the experimentally observed mag- 
netic structure. Spins are indicated by arrows, Mn d 4 sites 
by double-lobed orbitals, 2 ~ ions by shaded circles and 0~ 
ions by black circles. 



Mn-0~-Mn and Mn-0~Mn-0~Mn chains constitute 
the polarons for x = 1/2 and x = 2/3, respectively. Mn 
magnetic moments are in the range 3.8 to 3.9 (1b for 
either compound while oxygen ions within the polarons 
have charges closer to than 2 ~ and magnetic mo- 
ments of -0.7 /x b ■ All other oxygen ions are essentially 
2 ~ ions. Each polaron has a net magnetic moment of 
7.0 /ib, in good agreement with recently reported mag- 
netisation measurements from electron holography and 
Fresnel imaging 17 . 

UHF calculations predict that magnetic coupling 
within polarons is strong and ferromagnetic (FM) while 
coupling between polarons is much weaker and can be an- 
tifcrromagnetic (AF) or FM; these observations lead to 
a natural explanation for the observed CE or A-type AF 
magnetic ground states of manganites with x = 1/2 4,8 
and the magnetic ground state of Lao.33Cao.67Mn03, 
which contains polarons with the magnetic moment 
roughly aligned along the polaron axis 6 . Given that 
the magnetic coupling within the polarons is strong and 
FM, the appropriate model Hamiltonian for these sys- 
tems at low temperature is a Heisenberg Hamiltonian 
on a triangular lattice where each polaron is a single 
magnetic entity, rather than the conventional double ex- 
change Hamiltonian. 

Results of UHF our calculations are at odds with den- 



sity functional theory (DFT) calculations in the local spin 
density approximation (LSDA) reported recently 18,19 
and with LSDA calculations done as part of this work 
using the same electronic structure code 20 in that UHF 
calculations predict a magnetic moment on oxygen ions 
within the polarons whereas DFT calculations do not. 



II. UHF CALCULATIONS 

All electron UHF calculations were performed for 
Lao.sCao.sMnOs and Lao.33Cao.67Mn03 using the low 
temperature structures refined using x-ray data by 
Radaelli and coworkers 21 in the full (80 and 120 atom) 
magnetic unit cells. 

Spin densities of Lao.sCao.sMnOs (Fig. 3) and 
Lao.33Cao.67Mn03 (Fig. 4) in the ac planes of the crystal 
structures clearly show the polaronic nature of the elec- 
tronic structures, including the magnetic moment of the 
0~ ions opposed to those of the neighbouring Mn ions. 
The magnetic unit cells of either compound each contain 
four polarons in the ac planes shown. 




FIG. 3. UHF spin density plot for Lao.sCao.sMnOs in the 
CE-AF state. 



Magnetic moments and Mulliken populations on Mn 
and oxygen ions in the polarons are given in Table I. 

Electron holography and Fresnel imaging 17 measure- 
ments show that below the AF Neel temperature the x 
= 1/2 compound actually consists of both FM and AF 
domains and that the magnetic moment per Mn ion in 
the FM domains was 3.4 ± 0.2 /xg. UHF calculations 
predict a net magnetic moment of 7.0 \ib for the po- 
laron, i.e. the polarons are fully spin polarised, in agree- 
ment with these measurements of the magnetisation of 
FM domains 17 . Earlier neutron scattering data for this 
compound (Table I) indicated magnetic moments some- 
what smaller than these fully polarised values 5 but the 
trend, where larger moments arc found on the sites de- 
noted Mn 3+ , is reproduced. 
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FIG. 4. UHF spin density plot for Lao.33Ca .67Mn0 3 in AF 
state. 

TABLE I. Charges and magnetic moments in 
Lao.sCao.sMnOa from Mulliken populations and experiment. 





Mn moment /is 


O moment hb 


UHF (CE-AF) a 


3.78, 3.91 


-0.67 


Expt. ab 


2.57, 2.98 




Expt. c 


3.4 ± 0.2 






Mn charge 


charge 


UHF (CE-AF) a 


2.18, 2.17 


-1.24 



a Two values for the Mn moment are quoted as the Mn 
ions in the polaron are inequivalent 
b Neutron scattering, Ref. 5 

c Electron holography, value obtained for FM domain per 
Mn ion, Ref. 17 



Mulliken populations from UHF calculations on 
Lai-zCa^MnOs indicate total Mn ion populations in the 
range 2.25 (x = 0) to 2.13 (x = 1) with a monotonic vari- 
ation for intermediate values of x. Mn d populations have 
an even smaller relative variation across the range of x 
with a d population of 4.66 in LaMn03, a range of d pop- 
ulations from 4.72 to 4.73 in Lao^Cao.sMnOa, and 4.70 
in CaMnC>3. This is consistent with a bonding picture 
in which the t2 9 shell on each Mn ion is half-filled and a 
pair of c g orbitals is combined with a set of four empty 
sp 3 orbitals to form a set of equivalent d 2 sp 3 octahedral 
hybrid orbitals for polar-covalent Mn-0 bonding. The 
consistency of both the Mn ion population and d pop- 
ulation across the range of values of x contradicts the 
conventional double exchange picture where Mn ions are 
assumed to have their formal Mn 3+ and Mn 4+ charges. 

Total energies were calculated for the FM state and 
four different low energy AF states of the x = 1/2 com- 
pound. Low energy states are found by flipping entire 
polaron magnetic moments; flipping the spin of one of 
the Mn ions in a polaron results in an increase of the 
total energy by around 400 meV. Energies of the various 



low energy magnetic states of the x = 1/2 compound are 
well-fitted by an Ising-like, nearest-neighbour Hamilto- 
nian of the form given in Eq. 1 . 



H 



<ij> 



i ■ S- 



S 2 



(1) 



Labelling of exchange couplings is shown schematically 
in Fig. 5 and exchange constants obtained by fitting total 
energies of UHF calculations are given in Table II. The 
spin magnitude in Eq. 1 was chosen to be S = 2. 




FIG. 5. Schematic illustration of exchange couplings be- 
tween polarons in a Heisenberg Hamiltonian description. 



Exchange couplings within the plane shown in Fig. 5 
along zig-zag chains are FM and roughly equal in mag- 
nitude whereas they are AF perpendicular to the zig-zag 
chains. Exchange coupling between polarons in differ- 
ent planes is AF. The magnetic ground state is expected 
to be CE-type when Jafi exceeds (Jfmi + Jfmi)I^ in 
magnitude and A-type otherwise. Both are observed to 
be the magnetic ground state for various x = 1/2 com- 
pounds, depending on the A and B ion types 8 . 

TABLE II. Exchange coupling constants derived from 
Lao.sCao.sMnOs UHF calculations. 



Exchange constant (meV) 


Jafi 


5 


JaF'2^ 


8 


Jfmi 


-14 


JfM2 


-12 



a Average value over two interplanar couplings in crystal- 
lographic unit cell 

Our calculations actually predict an A-type AF ground 
state whereas the CE-AF structure is the ground state 
in Lao.sCao.sMnOa 10 . However, in similar UHF calcula- 
tions on LaMn03 22 , it was found that AF couplings are 
significantly underestimated compared to either experi- 
ment or results of configuration interaction (CI) cluster 
calculations, whereas FM couplings are in better agree- 
ment with both. 



III. DISCUSSION 

The polaronic picture for manganites with x > 1/2 
which has been presented is consistent with a wide range 
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of experimental data and allows various observations, 
such as the unusual magnetisation in Lao.33Cao.67Mn03, 
to be explained: it accounts for the observation of A- 
AF or CE-AF ground states for various combinations of 
counterion in Ao.sBo.sMnOa 4,8 ; it is consistent with full 
spin polarisation of Mn ions in Lao.sCao.sMnOs 17 . 

One must be cautious in using UHF calculations to es- 
timate magnitudes of magnetic moments on the oxygen 
ions. An analogy can be drawn between electron transfer 
from the O ion in a polaron to a neighbouring Mn ion and 
separation of the electron pair in a hydrogen molecule as 
the proton-proton distance is increased above the molec- 
ular equilibrium bond length: UHF wave functions for 
hydrogen molecules at large bond distances consist of a 
spin-up electron on one proton and an spin-down elec- 
tron on the other; the additional configuration in which 
spins on protons have been exchanged is omitted, owing 
to the simple functional form of the wave function. UHF 
(and LSDA) wave functions incorporate strong electron 
correlation effects at the cost of a proper treatment of 
electron spin. Spin symmetry is restored in CI methods, 
but these can only be applied to finite clusters owing to 
the large number of electronic configurations encountered 
in that approach. CI calculations on the pair of MnOg 
octahedra in the polaron in Lao.sCao.sMnOs show that 
the overall spin configuration is a linear combination of 
several determinants 23 with a magnetic moment on the 
central oxygen ion. The Lao.sCao.sMnOs polaron elec- 
tronic structure therefore resembles that in a'-NaV2 5 24 ; 
it is unlike the electronic structure which results in FM 
coupling in LaMnOs, which is very well described by a 
single spin configuration 22 . 

An obvious question which arises within the polaron 
picture is, 'What happens when electrons are added to 
the x = 1/2 phase?' That is, 'What is the electronic 
structure predicted to be in the FM region of the phase 
diagram with 0.2 < x < 0.5?' Since the bottom of the 
conduction band in Lao.sCao.sMnOs is comprised largely 
of vacant oxygen 2p states on the 0~ ion in the polaron, 
one would naturally expect to add the extra electron here. 
The spin density plots in Figs. 3 and 4 show that 00 in 
the x = 1/2 compounds is parallel to the polaron axis, 
whereas it is T-shaped in LaMnOs (long Mn-0 bonds 
in JT distorted octahedra in LaMnOs are perpendicu- 
lar to each other). If the extra electron is added to the 
polaronic oxygen site, it is expected that there will be 
reorientation of e g orbitals (Fig. lc) and adjustment of 
Mn-O-Mn bond lengths. 00 in the vicinity of the added 
electron is expected to resemble that in LaMnOs and one 
might expect to nucleate a small FM patch as further 
electrons are added. 
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